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Abstract 

In order to conduct a thermal analysis of heat transfer to liquid hydrogen near the 
critical point, an accurate understanding of the thermal transport properties is required. A 
review of the available literature on hydrogen transport properties identified a lack of 
useful equations to predict the thermal conductivity and viscosity of liquid hydrogen. 
The tables published by the National Bureau of Standards were used to perform a series 
of curve fits to generate the needed correlation equations. These equations give the 
thermal conductivity and viscosity of hydrogen below 100 K. They agree with the 
published NBS tables, with less than a 1.5 percent error for temperatures below 100 K 
and pressures from the triple point to 1000 KPa. These equations also capture the 
divergence in the thermal conductivity at the critical point. 
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Introduction 

Typical analytical forms of the equation of state predict finite values for the 
specific heats C p and C v , the thermal conductivity, and the speed of sound at the critical 
point. It is well established by experimental results that the specific heats and the thermal 
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conductivity diverge to infinity as the critical' point is approached. Garland and Williams 
[1] studied the speed of sound in xenon near its critical point. The speed of sound varies 


as } . Therefore, as C v goes to infinity near the critical point, the speed of sound will 

go to zero. They show that although the speed of sound approaches zero at the critical 

(T-T c ) _ 9 

temperature, it is still larger than 10 m/s for >10 . The speed of sound going to 

T 


zero at the critical point is common to all fluids. 


The use of empirical equations of state is required to adequately model the change 
in the fluid properties near the critical point. Hendricks et al. [2] developed a computer 
program to provide the thermophysical properties of ten different cryogenic fluids. The 
program did not attempt to model the anomalous behavior near the critical point, and 
therefore could not be utilized for research in critical point heat transfer. Emanuel [3] 
applied critical point theory to predict the thermodynamic properties of a real fluid near 
the critical point. The theory is applicable to any fluid and accurately predicts the 
divergence of the properties of interest at the critical point. The use of critical point 
theory needs further investigation and is beyond the scope of this research. Nwobi et al. 
[4] recently used molecular dynamics theory tosjiredict the thermodynamic properties of 
supercritical oxygen. The approach accurately predicts the divergence of certain 
important properties near the critical point, but is much too time intensive to be of any 
practical use in numerical simulations. 


In 1 980, the National Bureau of Standards (NBS) [5] created a computer program 
for the NASA Johnson Space Center to predict the thermophysical properties for fifteen 
cryogenic fluids. The program accurately reproduces the divergence of the fluid 
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properties at the critical point. The output from the program was compared to the 
published tables from the NBS with very close correlation, often within 1% of the 
tabulated values. This program did not include the thermal conductivity and viscosity of 
hydrogen. A review of previous reports from the NBS [6,7] on the thermal conductivity 
and viscosity of liquid hydrogen failed to provide any useful equations to predict these 
two transport properties. A decision was made to perform a series of curve fits using the 
published values from the NBS. 

Thermal Conductivity 

The NBS tables [8] were used to curve fit the thermal conductivity of gaseous and 
liquid hydrogen. The curve fits are valid for temperatures from the triple point to 100 K 
and pressures to 1000 KPa. For temperatures below 100 K, the thermal conductivity data 
were fitted to a modified Enskog theory using Eq. (1). 

k = exp(Z + exp(B) * p + C * p 2 ) (1) 

where the coefficients are given as 

Y = ln(T) 

Z = X, + X 2 * Y + X 3 * Y 2 + X 4 * Y 3 + X 5 * Y 4 

b = x 6 +x 7 *y+x 8 *y 2 +X 9 * y 3 +x io * y 4 +x n *y 5 
C = x i2 + X 13 /T + X u /T 2 + X 15 /T 3 + x |6 /T 4 

This form of the modified Enskog equation was used to predict the thermal 
conductivity of helium [9]. Since the properties of hydrogen are close to those of helium 
it was decided to use the same form of the modified Enskog equation to predict the 
thermal conductivity of hydrogen. The temperature, T, is in degrees Kelvin and the 
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density, p, is in g/cm 3 . Thermal conductivity, k, is given in W/m-K. The resulting curve 
fit coefficients are given as 


XI = 0.383546x10 + 01 
X3 = 0.611632x10 + 01 
X5 = 0.922390x10-01 
X7 =-0.682941x10 + 02 
X9 =-0.884876x10 + 0 
XI 1 = - 0.568872x10 -01 
X13 = - 0.300763x10 + 04 
X 1 5 = 0. 1 06809x 10 + 07 


X2 = - 0.123663x10 + 02 
X4 =-0.123492x10 + 01 
X6= 0.569497x10 + 02 
X8= 0.348056x10 + 02 
X10= 0.112028x10 + 01 
X12= 0.473430x10 + 02 
X14 = 0.776719x10 + 05 
X16= 0.888314x10 + 06 


The average error was 0.046% and the standard deviation was 4.35%. A three- 
sigma spread on the data would give an estimated 1 3% error. The maximum error differs 
by 99% from the published tables near the critical point. Far away from the critical point 
the error is generally less than 0.5%. The results from Eq. (1) showed very good 
correlation at the higher densities. 
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The difference between the values predicted by Eq. (1) and the published values 
demonstrated a seventh order polynomial relationship based on density for constant 
temperature conditions. This relationship leads to the use of seventh order polynomials 
for density and cubic polynomials for temperature to correct for the lower density 
conditions. The resulting form of the correction term is given by Eq. (2). 


Corr = A + B*p + C*p 2 +D*p 3 +E V + F*p 5 +G*p 6 + H V (2) 


where the coefficients are given as 

A = X, + X 2 *T + X 3 *T 2 +X 4 *T 3 
B = X 5 +X 6 * T + X 7 * T 2 + X g * T 3 
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C = X 9 + X l0 * T +X n * T 2 + X 12 *T 3 


D = X| 3 + X 14 * T + X 15 * T 2 + X l6 * T 3 

e = x 17 +x, 8 *t + x, 9 *t 2 + X 20 *t 3 

F = X 21 + X 22 * T + X 23 * T 2 + X 24 * T 3 

G = X 25 + X 26 * T + X 27 * T 2 + X 28 * T 3 

H = X 29 + X 30 *T + X 31 *T 2 + X 32 *T 3 

The coefficients of the density correction for different temperature and density 
ranges are given in Table 1. The final form of the equation for thermal conductivity is 
the sum of Eqs. (1) and (2). The results are summarized in Table 2. 

The critical point enhancement was applied to the thermal conductivity in the 
temperature range 32.95K <T <33K and the density range 0.02 <p< 0.04g/cm 3 . The 
critical point enhancement was calculated by Eq. (3). 

Crit = 41 .20535 * exp(-6675067* (p - 0.03 1428) 2 - 2.0xl0 7 * (T - 3? 976) 2 ) (3) 

Figure 1 shows the maximum error between the calculated values and the values 
given by the NBS tables for temperatures below 100 K. The data fall below a 1.6% error, 
with the majority well below 1%. The curve fit was then compared to the tabular values 
for three significant pressures: atmospheric, critical, and super-critical. Figure 2 
compares the predicted thermal conductivity to the published values for hydrogen at one 
atmosphere pressure. The predicted values closely match the published values for both 
the liquid and gaseous phases. The correlation at the critical pressure is presented in 
Figure 3. The curve fit results accurately capture the divergence at the critical point. The 
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divergence is limited to a very narrow temperature range. Finally, the correlation at a 
supercritical pressure is shown in Figure 4. Here the pressure was set to 1500 KPA or 
1.16 times the critical pressure. From the shape of the curve it is obvious why three 
density ranges were required to provide an acceptable curve fit. The resulting equations 
still capture the divergence near the critical point. The magnitude of the critical point 
divergence decreases as the pressure is increased. 

Viscosity 

The same NBS tables [8] were used to curve fit the viscosity of hydrogen. The 
curve fits are valid for temperatures from the triple point to 1 00 K and pressures to 1 000 
KPa. 

In the temperature range 65K < T ^ 100K the data were fitted using Eq. (4). 

p = A+B*p+C*p 2 +D*p 3 +E*p 4 +F*p 5 (4) 

where the coefficients are given as 

A = 0. 1 1 34 1 4x 1 O’ 06 + 0.524985x1 0" 07 * T - 0. 1 1 5 1 24x1 O' 09 * T 2 
B =0.232625x1 O' 03 - 0.453684x1 O' 05 * T + 0.244673x1 O' 07 *T 2 
C = -0.245986x1 0" 01 + 0.551425xl0~° 3 * T - 0.295357xl0"° 5 * T 2 
D = 0.10104 lxl0 +O1 - 0.218787x10“°' * T + 0.1 15788x1 O' 03 * T 2 
E = -0.1 56836x1 O' 102 + 0.338478xl0 + °° *T - 0.1771 12x1 0" 02 *T 2 
F = 0.856032x1 0 + ° 2 -0.180672xl0 + °‘ *T + 0.933 153x1 0"° 2 *T 2 
The temperature, T, is in degrees Kelvin and the density, p, is in g/cm 3 . 
Viscosity, p, is in Poise. The average error was 0.0013% and the standard deviation was 
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0.32%. A three-sigma spread on the data would give an estimated 0.96% error. The 
maximum error for this range was 1 .62%. 

In the temperature range 46 K < T < 65 K the data were fitted using Eq. (5). 

p = A+B*p+C*p 2 +D*p 3 +E*p 4 +F*p 5 + G*p 6 +H*p 7 (5) 

where the coefficients are given as 

A = 0.542644x1 0" 06 + 0.3 84834x1 O' 07 * T + 0.373781x10"" * T 2 

B = -0.382474x1 0" 03 + 0. 1 56988x1 O' 04 *T-0.1 44832x1 O' 06 * T 2 

C = 0.648002x10'°' - 0.262206x1 O^ 2 * T + 0.249320x1 0'° 4 * T 2 

D = -0.343722xl0 + °' + 0. 1521 63x1 0 + °° *T-0. 149907x1 O' 02 * T 2 

E = 0.82431 2x1 0 +02 - 0.404907x1 0 +O1 * T + 0.415405x10"°' * T 2 

F = -0.9061 90x1 0 + ° 3 + 0.528978x1 0 + ° 2 * T - 0.573267x10^° * T 2 

G = 0.376977x1 0 + ° 4 - 0.319712xl0 +O3 * T + 0.377532xl0 + °' * T 2 

H = 0 467651x1 0 + °° +0.656682x1 0 + ° 3 * T - 0.9 1 6660x 1 0 + °' * T 2 

The average error was 0.031% and the standard deviation was 0.25%. A three- 
sigma spread on the data would give an estimated 0?78% error. The maximum error for 
this range was 1.53%. 

For T <46K and p < 0.075g/cm 3 the data were fitted using Eq. (6). 

p LD = A+B*p+C*p 2 +D *p 3 +E*p 4 (6) 

where the coefficients are given as 

A = X, + X 2 *T + X 3 *T 2 +X 4 *T 3 +x 5 *t 4 
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b = x 6 + x 7 *t + x 8 *t 2 +x 9 *t 3 +x I0 *t 4 


C = X„ + X I2 *T + X, 3 *T 
D = X, 6 + X 17 * T + X l8 * T 
E = X 21 +X 22 *T + X 23 * 

where, 

XI = 0.189281x10-06 
X3 = 0.277990x10-08 
X5 = 0.599499x10-12 
X7 =-0.219076x10-04 
X9 =-0.255997x10-08 
XI 1= 0.732495x10-01 
X13= 0.302109x10-03 
XI 5= 0.380544x1007 
X17= 0.266919x10 + 00 
X19= 0.187012x10-03 
X21= 0.229229x10 + 02 
X23= 0.891753x10-01 
X25= 0.102719x10-04 


+ X I4 *T 3 + X I5 *T 4 
+ X 19 *T 3 +X 20 *T 4 
+ X 24 * T 3 + X 25 * T 4 


X2 = 0.130968x10-07 
X4 =-0.721669x10-10 
X6 = 0.266869x10-03 
X8 = 0.541725x10-06 
X10 = - 0.3 17459x10- 10 
X12 = - 0.744391x10 -02 
X14 = - 0.5533 16x10 -05 
XI 6 = -0.262958x10 + 01 
XI 8= -0.1 05438x10 -01 
X20=- 0.1 24540x10 -05 
X22 = - 0.228625x10 + 01 
X24 = - 0.1 561 64x10 -02 


The average error was 0.078% and the standard deviation was 0.54%. A three- 
sigma spread on the data would give an estimated 1.70% error. Thejjnaximum error for 
this range was 2.39%. 


This same equation was used to extrapolate into the higher density region. A 
density correction term was calculated to reduce the magnitude of the error. The final 
form of the equation for viscosity is the sum of the two terms, p = p LD + Phd (Pm>i s 
obtained by substituting the correction coefficients into Eq. (6)). The coefficients of the 
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density correction term are as given below for different temperature and density ranges. 
The effects of correction for these cases are given in Table 3. 

In the temperature range 40K <T< 46K and p>0.075g/cm\ the definition of 
the coefficients was changed to reflect a linear relationship in temperature in place of the 
previous quadratic temperature relationship. The coefficients are now defined as 

A = 0.296606x1 0" 01 - 0.599524x1 0^ 3 * T 

B = -0.14 1972x10 401 + 0.285903x10^' *T 

C = 0.254884x1 0 +02 -0.511 198xl0 +0 ° * T 

D = -0.20352 1x10 +03 + 0.406355xl0 + °' * T 

E = 0.61 0254x1 0 +03 - 0.121242xl0^ 2 * T 


36K<T<40K and p>0.075g/cm 3 

A = 0.538007x10"°' - 0.1 20407x1 O' 02 * T 
B = -0.253 839xlO +01 + 0.566168x10^' * T 
"3 C = 0.4491 75x1 0 +02 - 0.9981 52 xlO +0 ° * T 

D = -0.353461xl0 +03 + 0.72285xl0 +O1 * T 
E = 0.1 04424x1 0 +04 - 0.230085x10^ * T 


33K <T <36K and p>0.075g/cm 3 

A = -0.275020x10'°' +0.1 102 12x1 0^ 2 * T 
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B =0.1 17525xlO +01 -0.487800x1 O' 01 *T 
C = -0.1 843 17x10 402 +0. 800791 xlO +0 ° *T 
D = 0. 124700x1 0 +03 - 0.576384X10 401 * T 
E = -0.302676x1 0 403 + 0.1 52884x1 0 +02 * T 


22K <T < 33 K and p > 0.075 g/ cm 3 

A = 0.282 191x1 0 -01 - 0.155902xl0"° 2 * T + 0.322665x1 0” 04 * T 2 
B = -0.131286x1 0 +01 + 0.685339x10-°' * T - 0.13991 lxl O -02 * T 2 
C = 0.231542xl0 + ° 2 - 0.1 1 3302x1 0 +O1 * T + 0.226347x1 O' 01 * T 2 
D = -0. 1 84 1 77x 1 0 +03 + 0.83944 lxlO +0 ‘ * T-O.l 62347x1 0 +0 ° * T 2 
E = 0.559566x1 0 +03 -0.236825xl0 +O2 * T + 0.437404x1 0 400 *T 2 


20 K < T < 22K and p > 0.077 g / cm 3 

A = -0.1 9653 8x1 O' 01 + 0.1 288 18x1 O’ 02 * T 
B = 0.90154Oxl0 +o ° - 0.6 1 3574x1 0 -01 * T 
C = -0. 1521 79x1 0 +02 + 0.108723xl0 +O1 * T 
D = 0.1 1 1 263x1 0 +O3 - 0.847962x1 0 + °‘ * T 
E = -0.293960x1 0 + ° 3 + 0.245057xl0 + ° 2 * T 

17K < T< 20K 
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A = -0.446732x1 0"°' + 0.252585xl0 r ° 2 * T 
B = 0.224243xl0 +t)1 - 0.127828xl0 + °° * T 
C = -0.420352x1 0 +02 + 0.24 1908x10 401 * T 
D = 0.348548x1 0 403 - 0.202838x1 0 +02 * T 
E = -0. 1 07787x 1 0 +04 + 0.656 1 7x 1 0 +02 * T 

T < 17K and p > 0.075g/cm 3 

A = -0.853063x10'°' + 0.976589xl0'° 2 * T - 0.247566x1 O' 03 * T 2 
B =0. 321564 xlO +01 -0.35681 5x1 0 + °° * T + 0.813 185xlO’ 02 * T 2 
C = -0.392538xl0 4 ° 2 + 0.407 149xlO +0 ' * T- 0.682502x1 O' 01 *T 2 
D = 0.1 39394x1 0 + ° 3 -0.1 07 156x1 O 402 * T - 0.174518x10^° * T 2 
E = 0.1 94 107x10 403 - 0.438586x1 0 + ° 2 *T + 0.301042xl0 + °' *T 2 

Figure 6 shows the maximum error between the calculated values and the NBS 
table values for viscosity at temperatures below 1 00 K. The data fall below a 5% error 
with the majority below 2%. The curve fit was then compared to the tabular values for 
three significant pressures: atmospheric, critical, and supercritical. Figure 7 compares the 
predicted viscosity to the published values for hydrogen at one atmosphere pressure. The 
predicted values closely match the published values for both the liquid and gaseous 
phases. The correlation at the critical pressure is presented in Figure 8. Viscosity does 
not show a divergence at the critical point. Finally, the correlation at a supercritical 
pressure is shown in Figure 9. Here, the pressure was set to 1500 KPA or 1.16 times the 
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critical pressure. From the shape of the curves it is obvious why only a low density and a 
high-density region were required to provide an acceptable curve fit. 


Conclusion 

The equations described in this paper provide an easy and accurate method of 
predicting the thermal conductivity and viscosity of hydrogen for temperatures below 1 00 
K and pressures up to 1000 KPa. The FORTRAN program created by the NBS [5] was 
modified to include the equations described above. The computer program was then 
employed to perform a comparison of the pressure response of liquid hydrogen to the 
pressure response of liquid oxygen at their respective critical points. 
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Table 1 - Thermal Conductivity Curve Fit Coefficients 



40K < T < 75K 

33K < T < 40K 

32K < T < 40K 

T < 32K 

T < 22K 


Curve fit 1 

Curve fit 2 

Curve fit 3 

Curve fit 4*3 

Curve fit 5 

X, 

0.359343xl0' 2 

0.610708x10'' 

0.996513xl0" 2 

0.231879x10'' 

-0.302237x10"' 

X 2 

-0.320203x1 0' 3 

-0.480935xl0' 2 

-0.670513x10"' 

-0.29 1474x1 0' 2 

-0.654848x10"° 

X 3 

0.684352x1 0' 5 

0.122292x1 O' 3 

0.186958x10"° 

0.1 13756x1 0' 3 

0.703479x10'' 

X4 

-0.429492x1 (T 7- 

-0.101920xl0’ 5 

-0.188593x1 O' 2 

-0.1 44868x1 0' 5 

-0.1 45 155x1 O' 2 

X 5 

-0.657379x1 0 +0 

-0.904625x1 0" 2 

-0.283875x1 0 +4 

-0.571286x10"' 

0.169052xl0 +3 

x 6 

0.772415x10'' 

0.688178x10"' 

0.1 13848xl0" 3 

0.882118x10"° 

-0.291281x10"' 

X 7 

-0.196103xl0' 2 

-0.170487x10"° 

-0.448288x10"' 

-0.401339x10'' 

-0.474153x10"° 

X 8 

0.142554X10- 4 

0.137255xl0' 2 

0.646446x10'' 

0.572794x1 0' 3 

0.403391x1 O' 2 

X 9 

0.147837xl0" 3 

0.21 7502x1 0" 5 

0.149707xl0 +5 

-0.556360x1 0" 3 

-0.129162xl0 +3 

X,o 

-0.180362xl0" 2 

-0.155565xl0" 4 

0.384486xl0" 4 

0.71761 lxl0 +1 

0.426 176xl0 +2 

X,, 

0.490849x1 0 +0 

0.350012xl0" 2 

-0.423525xl0" 2 

0.215320xl0 +1 

-0.556174xl0 +1 

X,2 

-0.376166xl0' 2 

-0.241006x10"° 

-0.403244x1 0 + ° 

-0.584037x10'' 

0.353707x10"° 
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X , 3 

0.291544x10"® 

-0.105106xl0 +7 

-0.11 0272x1 0" 7 

0.524678x10"° 

-0.179205xl0' 3 

X , 4 

-0.506066x1 0" 3 

0.570871x10"® 

-0.525837xl0" 5 

-0.434403x10"® 

0.185432xl0 +4 

X, 5 

-0.902909x1 O’’’ 1 

-0.514534xl0" 3 

0.719500xl0" 3 

0.101659xl0" 4 

-0.689358x1 O' 2 

X , 6 

0.148432x10"° 

-0.651237xl0 +1 

-0.487915x10"' 

-0.495007x10"' 

0.537557x10"° 

X , 7 

-0.120000xl0" 7 

0.129222xl0 +2 

0.391219x10"® 

-0.668844x10"' 

-0.995309x10"® 

X , 8 

0.241617x10"® 

0.122469xl0 +7 

-0.710429x10"° 

0.224704x1 0" 6 

0.974862x1 0" 4 

X, 9 

0.272422x1 0" 3 

-0.861367xl0 + ® 

0.1 1 1 144xl0" 5 

-0.346450x1 0" 4 

-0.305666xl0 +3 

X 20 

-0.542153x10"' 

0.136533xl0 +4 

0.251 145xl0" 3 

-0.150762xl0" 3 

-0.453023x1 0" 2 

X 21 

0.1 13418x10"® 

0.697990x10"' 

-0.324325x1 0" 9 

-0.41 60 1 Ox 1 0 +1 

-0.512421xl0" 7 

X 22 

0.257187x10"® 

-0.644070x1 0" 2 

0.315749xl0" 7 

0.314617x10"' 

0.210300xl0 + ° 

X 23 

-0.1 23460x1 0" 5 

0.531457xl0 +6 

0.107766x10"° 

-0.279761x10"° 

0.113568x10"® 

X 24 

0.141 239x1 0 +3 

-0.933265xl0 +4 

-0.745936x1 0" 4 

0.162868xl0" 4 

0.416300xl0" 3 

X 25 

0.504982x10"° 

0.000000x10"° 

0.1 37662x1 0 +3 

0.000000x1 o +0 

0.765159xl0 + ® 

X 26 

-0.357369xl0" 7 

-0.20421 Ox 10 +2 

0.890544x10"® 

-0.402382x10"' 

-0.266769x1 0 +7 

X27 

0.146900xl0" 6 

-0.470943x1 0 +2 

-0.287183xl0" 7 

-0.117582x10"' 

-0.146626x10"° 

X 28 

-0.141778xl0 +4 

-0.197265x10"° 

0.655508x10"® 

-0.9041 llxlO +l 

-0.144348xl0 +4 

X29 

-0.489403x1 0" 2 

0.000000x1 o +0 

0.668644xl0" 2 

0.000000x10"° 

-0.379681xl0" 3 

X30 

0.162709xl0" 2 

0.000000x1 o +0 

-0.340703x1 0" 3 

0.105514xl0’ 7 

-0.3 1972 lxl 0 + ® 

X 31 

-0.1 5 1925x1 0" 6 

-0.553572xl0" 2 

-0.540557xl0" 7 

0.207713xl0 + ° 

0.234894x1 0 +7 

X32 

0.245402x1 0" 4 

0.3 1442 lxlO" 7 

0.120796xl0 +4 

-0.5541 18xlO +1 

-0.231774xl0 + ® 


-3 


14 


Table 2. The effect of density correction for thermal conductivity 


Case 

DC 

AEBC 

AEAC 

STBC 

SDAC 

TSBC 

TSAC 

MEBC 

MEAC 

1 

No 

0.097 

- 

0.40% 

- 

1.3% 

- 

0.88 

- 

2 

Yes 

0.13 

0.0059 

1.74 

0.25 

5.35 

0.75 

11.56 

1.6 

3 

No 

0.22 

- 

0.25 

- 

0.97 

- 

0.88 

- 

4 

- 

0.079 

- 

0.27 

- 

0.89 

- 

1.60 

- 

5 

Yes 

0.083 

0.013 

8.02 

0.16 

24.24 

0.49 

80.35 

0.64 

6 

Yes 

3.63 

0.72 

10.31 

0.12 

34.56 

1.08 

51.98 

0.88 

7 

No 

0.34 

- 

0.20 

- 

0.94 

- 

0.92 

- 

8 

Crit 

97.2 

0.098 

- 

- 

- 

- 

- 

- 

9 

- 

0.12 

- 

0.35 

- 

1.17 

- 

0.92 

- 

10 

Yes 

0.98 

0.0002 

6.03 

0.19 

19.07 

0.57 

28.7 

0.93 

11 

No 

0.022 

- 

0.26 

- 

0.80 

- 

0.72 

- 

12 

Yes 

0.16 

0.0057 

1.09 

0.12 

3.43 

0.42 

2.47 

0.40 

13 

- 

0.019 

- 

0.19 

- 

0.59 

- 

0.93 

- 


1- 75<T<100 

2- 40<T<75 and p<0.06, Curve fit 1 

3- 40<T<75 and p>0.06 

4- 40<T<75 for all densities after correction 

5- 32<40<40 and p<0.03 143, not including the critical point. Curve fit 2 

6- 32<T<40 and 0.03 143<p<0.06, not including the critical point. Curve fit 3 

7- 32<T<40 and p>0.06 

8- Critical point 

9- 32<T<40 for all densities after correction 

10- T<32 and p<0.03143, Curve fit 4 
1 1 - 22<T<32 and p>0.06 

12- T<22 and p>0.06, Curve fit 5 

13- T<32 for all densities 
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Table 3. The effect of density correction for viscosity 


Case 

AEBC 

AEAC 

SDBC 

SDAC 

TSBC 

TSAC 

MEBC 

MEAC 

1 

12.79 

0.24 

11.47 

0.29 

47.2 

1.11 

43.65 

0.93 

2 

18.07 

0.46 

14.59 

0.49 

61.84 

1.93 

53.63 

1.72 

3 

20.99 

0.12 

16.7 

0.80 

71.09 

2.52 

58.51 

4.36 

4 

18.28 

0.28 

14.22 

0.64 

60.94 

2.20 

56.05 

5.04 

5 

12.98 

2.17 

8.37 

0.47 

38.09 

3.58 

31.07 

3.06 

6 

7.33 

0.90 

5.62 

0.27 

24.19 

1.71 

20.27 

1.57 

7 

3.16 

0.50 

1.96 

0.11 

9.04 

1.07 

9.60 

0.74 


1 - 40<T<46 and p>0 .075 

2- 36<T<40 and p>0.075 

3- 33<T<36 and p>0.075 

4- 22<T<33 and p>0.075 

5- 20<T<22 and p>0.077 

6- 17<T<20 and p>0.076 

7- T<17 and p>0.075 


DC: Density Correction, % 

AEBF: Average Error Before Correction, % 
AEAC: Average Error After Correction, % 
SDBC: Standard Deviation Before Correction, % 
SDAC: Standard Deviation After Correction, % 
TSBC: Three-sigma Before Correction, % 
TSAC: Three-sigma After Correction, % 

MEBC: Maximum Error Before Correction, % 
MEAC: Maximum Error After Correction, % 
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Error, % 


Fig. 1 Maximum thermal conductivity error for 100 K and below 
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Thermal Conductivity, W/m-K 


Fig. 2 Thermal conductivity at atmospheric pressure (101.325 KPa) 
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Fig. 3 Thermal conductivity at the critical pressure (1292.8 KPa) 
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Thermal Conductivity, W/m-K 


Fig. 4 Thermal conductivity at the supercritical pressure (1500 KPa) 
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Error, % 


Fig. 5 Maximum viscosity error for 100 K and below 
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Viscosity, millipoise 
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Fig. 6 Viscosity at atmospheric pressure (101.325 KPa) 
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Viscosity, millipoise 


Fig. 7 Viscosity at the critical pressure (1 292.8 KPa) 
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Viscosity, millipoise 


Fig. 8 Viscosity at the supercritical pressure (1500 KPa) 
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